Significant primordial star formation at redshifts z ~ 3 — 4. 
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Four recent observational results have challenged our understanding of high redshift galax- 
ies, as they require the presence of far more ultraviolet photons than should be emitted 
by normal stellar populations. First, there is significant ultraviolet emission from Ly- 
man Break Galaxies (LBGs) at wavelenghts shorter than 912A|T]. Second, there is strong 
Lyman a emission from extended "blobs" with little or no associated apparent ionizing 
continuum[2 . Third, there is a population of galaxies with unusually strong Lyman a emis- 
sion lines[3j. And fourth, there is a strong Hell (1640 A) emission line in a composite of 
LBGs[4j. The proposed explanations for the first three observations are internally incon- 
sistent, and the fourth puzzle has remained hitherto unexplained. Here we show that all 
four problems are resolved simultaneously if 10 — 30 percent of the stars in many galaxies 
at zfs3-4 are mainly primordial unenriched by elements heavier than helium ('metals'). 
Most models of hierarchical galaxy formation assume efficient intra galactic metal mixing, 
and therefore do not predict metal-free star formation at redshifts significantly below z ~ 5 
[5l l6l |7[ l8l l9j . Our results imply that micro mixing of metals within galaxies is inefficient 
on a ~ Gyr time scale, a conclusion that can be verified with higher resolution simulations, 
and future observations of the Hell emission line. 
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The continuum emission from LBGs can be explained by adding the continuous formation of massive stars 
to the standard initial mass function (IMF). In the absence of such a bias in the IMF, the continuum flux 
emerging just below < 912 A will be dominated by low-mass stars and will be significantly fainter than 
observed. The other three puzzles require a substantial increase in the total number of energetic ionizing 
photons, and we find they cannot be resolved by a bias of the IMF alone. To increase the number of 
ultraviolet photons further the temperature in the atmospheres of the stars must be increased. This is 
the case for metal-free (often called population-Ill) stars, which have significantly reduced opacity owing 
to the lack of CNO elements HOI ITT]. 

Lyman-continuum flux has been detected £Q in the composite spectrum of 29 Lyman Break Galaxies 
(LBGs) at a mean redshift of z sw 3.4. Although the authors of ref. 1 warn that sky background 
subtraction at the low observed flux level is a concern, when taken at face value, their result implies that 
a surprisingly large fraction of ionizing photons can typically escape from these galaxies. These authors 
found / csc > 0.5, where / csc is defined as the ratio of 900 A photons that escape the galaxy and the 1500 A 
photons that escape (the wavelengths are quoted in the rest-frame of the LBGs). Using more detailed 
fits of the observed spectrum, |12j and |13| found that the spectrum appears surprisingly blue, even if 
no additional Lyman continuum absorption is assumed to take place in the galaxy (beyond the dust 
absorption that also affects the spectrum at 1500 A). 

We compare the stacked composite spectrum of LBGs in JTJ with model galactic spectra, including metal- 
free stars. The intrinsic emission of the galaxies is modeled as the sum of "normal" star-formation (stellar 
populations with metallicity ratio > 10 _3 Z/Zq |14II15| L and metal-free, massive pop-Ill star-formation. 
For the "normal" component of the intrinsic stellar emission, we use the spectral synthesis models of |16| 
for a Salpeter IMF (with a slope a = 2.35) between 1 — 100M Q and a metallicity of Z = 0.4Z Q , and 
compute the spectrum for continuous star formation for 10 8 yr, which is the age for which a steady- 
state has been reached. For the 'metal-free' component, we use the models of ^7], and adopt the same 
Salpeter IMF, but limited to massive stars in the 100 — 500M Q range and with Z = 0. Note 
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that massive stars with > 100M Q have similar spectra, and our results should not depend significantly 
on the distribution of pop-Ill stellar masses above this threshold value. To simulate the observed, 
processed composite spectrum, we modify the emitted stellar template spectrum by including the effects 
of dust absorption, as well as the opacity of the intervening intergalactic medium (see Supplementary 
information) . 

As the top panel in Fig. ^ shows, in the 880A< A < 910A range, the mean observed flux is quite high, 
Fv/Fisoo — 0.064 ± 0.013; while the normal-only star formation model predicts a much lower value of 
Fu/Fxwa — 0.023 ± 0.01. Taken at face value, the data reveals an ionizing flux about three times higher 
than this model predicts, even if we do not include any continuum absorption by neutral gas in the 
galaxy. The data and the model do not agree at the > 3ct level. The model with pop-Ill star-formation 
with /popin = 0.5 almost generates the observed ionizing flux in the 880A< A < 910A range; producing 
a mean F v /F 150 o = 0.048 ± 0.02. 

Our conclusions regarding the LBGs apply to the bluest sub-sample. This is further supported by the 
analysis of the CIV absorption line associated with the stellar winds in massive stars that is observed in 
the spectra of LBGs. According to our prediction, the bluest LBGs should have weaker CIV absorption 
than the reddest LBGs, because pop-III stars are metal free and therefore contribute significantly fewer 
metals. This trend is indeed what is observed (see table 3 and section 5.3 in 0]). Detailed inspection of 
the numbers in table 3 of 0] shows that the weaker strength of the CIV line requires an abundance of 
~ 50% pop-III in the bluest LBGs. We note that the CIV profile in LBGs is a complex superposition of 
interstellar absorption, stellar-wind emission and absorption, and possibly nebular emission. Our models 
predict that the trend of decreasing CIV strengh in the bluest galaxies is present in the narrow interstellar 
absorption and any broad stellar component (which measures the total cumulative C output of all stars), 
but not necessarily in any of the other components. The fact that metal lines are observed in all LBGs 
indicates that LBGs are not composed purely of pop-III stars, but that there is a mixture in each galaxy 
of pop-III and normal stars. 
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There have been recent observations of a population of 35 Lya blobs at z = 3.1 in narrow and broad 
band spectra[2]- For about one third of the blobs, the amount of photons provided by a normal stellar 
population that fits the observed spectrum at wavelengths longer than 912 Ais not capable of producing 
the required number of ionizing photons to explain the observed Lya flux. If we assume case B (two 
Lya photons are produced for every three ionising photons from the stars) then a stellar population with 
Z = 10 _2 Z Q would produce at least 35% less photons than those necessary to account for the observed 
Lya flux. To reproduce the observed flux, about 30% of the stars would need to be metal-free. 

Figure [21 shows that about 20 — 30% of the blobs also have continuum colors consistent with the stars 
producing the ionizing Lya photons being metal-free. The fact that the stellar populations powering 
the extended emission are metal-poor is consistent with the fact that these sources are extended, and 
are presumably more likely to have been caught in the process of their assembly, when they are less 
metal-enriched overall . They could plausibly be identified, therefore, as the progenitors of LBGs. The 
rest-frame equivalent width for Lya ranges from 20 to 350 A. For models with constant star formation 
and Z > 10~ 5 Z Q the EW are 70 - 100 A. About one third of the blobs have EW larger than 100 A 
and therefore require metal-free stars. This is in agreement with what is found for the colors, although 
we note that there is little overlap among the individiual blue blobs and those with Lya excess. The 
good fit to the B — I colors, without including (case-B) reprocessing of ionizing radiation into Lyman a, 
implies that / esc is of the order of 100% because a smaller escape fraction would imply more flux on B 
and therefore significantly bluer B — I colors than those observed. The large escape fraction of ionizing 
radiation would be consistent with the value observed in the composite of LBG spectra. 

Strong Hell (1640 A) emission line from a composite of LBGs has been observed @J. They report 
equivalent widths for a sample of several hundred LBGs of about 2 A. This line is also very broad, with 
a FWHM of about 1500 km s _1 . Although this cannot be explained by stars without winds, it has been 
shown (e.g. |201 121| 1 that, aided by rotation, metal-free stars can loose significant amounts of mass 
and drive strong winds, similar to normal-metalicity stars. In this case, the He-ionizing photons from 
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the massive popIII stars will be reprocessed in the optically thick stellar outflows surrounding the star, 
making the line broad As discussed in 0|, normal WR stars cannot explain the Hell line observed in 
the composite LBG sample, because the requisite number of WR stars will overproduce the stellar CIV 
emission. 

However, metal-free massive stars with winds do not have this problem. As shown in l 223 there is a 
clear anti-correlation between the strength of the Hell 1640 line and the CIV line (see Fig. 13 in [22] for 
the case Z = 0.01 and lO -4 ^). For metal- free stars the equivalent width of the Hell 1640 line for a 
continuous star forming (1 M Q yr^ 1 ) model of metal-free stars at age 0.1 Gyr is 30 A, significantly larger 
than observed. So for this sample, the required amount of metal-free stars is only 10%. However, the 
equivalent width of the Hell 1640 line is very sensitive to the metallicity. If we use Z = 10 _7 Zq (instead 
of Z = 0), then the 912 A decrement doesn't change, so our conclusions above about the UV radiation 
from LBGs and the Lya blobs remain unchanged. But the He-ionizing flux, and the Hcl640 equivalent 
width is reduced by about a factor of 6 — 7. This implies a population of ~ 30% metal-free stars. 

Recent non-detections of the Hell 1640 line in a composite spectrum of 17 galaxies at z — 4.5 [2H], as 
well as from an individual Lya emitter at z = 6.33 |24) have provided tight upper limits, but still allow 
for a significant mass fraction (~ 50%) in Pop III stars. 

The Large Area Lyman Alpha survey 3. reports large Lya equivalent widths for a sample of 150 emitters 
at z » 4.5. The median value of the distribution is ~ 450 A and 60% have equivalent widths above 
240 A. No continuum emission was found in 30% of the emitters. However, the equivalent widths are 
spread over a wide range: from 10 to 5000 A. This implies that not a single type of stellar population 
can account the large spread of equivalent widths. The first suggestion that this large equivalent widths 
can be produced by very low metallicity stars is found in [3UJ. We find that for metal- free stars and 
constant star formation, the equivalent width is 1000 A for the IMF used in this paper and declines to 
300A for a Salpeter IMF ^7]. We performed a more detailed fit to the distribution in Using models 
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for metallicities Z > 10 _5 Z Q we find that bursts of ages greater than 10 7 years have equivalent widths 
as large as 300A whereas constant star formation models have about 200A. Therefore "normal" stellar 
populations can account for as much as 40% of the Lya emitters, while the rest need to be accounted for 
by metal-free stellar populations. 

Recently two stars with metallicity Z < 10~ 5 Z Q have been found in the Milky Way [21]. Although such 
low-mass stars ( ^ 0.8 M©) are not included in our model, they reveal the details of the earliest stages 
of chemical enrichment, and suggest an initial overproduction of Carbon. For stars with [Fe/H] < —3.5, 
40% have [C/Fe] > 1.0 and the trend increases to 100% for stars with [Fe/H] < -5. This trend can 
naturally arise in our model |20l l2~T] . 

Our results imply that micro-mixing of metals within galaxies is inefficient on a ~ Gyr time-scale, a 
conclusion that may be surprising, as current models of hierarchical galaxy formation assume efficient 
intra-galactic metal mixing, and therefore do not predict metal-free star formation significantly below 
z ~ 5 [S1IH1E1IH1E1- However, the presence of significant pockets of metal-free star-formation explains 
four puzzling observational results, and the hypothesis does not violate any other existing observations. 
The presence of metal-free gas pockets can be verified in future high resolution numerical simulations, 
and by future observations. Our models predict that in some of the LBGs there should be a detectable 
continuum shortward of 912A and a detectable individual Hell line. Long-integration observations of 
these individual galaxies in the future could be used tocharacterize their spectra, and to confirm the 
presence of metal- free stellar populations. In particular we predict that there should be a clear anti- 
correlation between Hel640 and CIV line strength and this should be a reliable test of the amount of 
metal- free stars at z ~ 3 — 5. 
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Figure 1: Lyman break galaxies stacked observed spectrum and best fit population synthesis 
models. The dotted curve corresponds to continuous normal star formation for 10 s years, and the 
dashed curve corresponds to a mix model where 50% of the stars (by mass) are massive, metal-free 
stars in the 100 — 500 M Q range. The solid curves show the composite spectrum in pQ. All fluxes are 
normalized to the flux at the emitted wavelength of 1500A. For reference, in the bottom panel, we show 
the optical depths assumed for the dust and the IGM. In the top panel, the horizontal dashed line shows 
the mean observed flux in the interval 880A< A < 910A. The inference is that the spectrum requires a 
significant contribution by metal-free star-formation, with a mass-fraction / p0 piii ~ 0.5. As discussed 
in |12j , changing the assumed metallicity of the normal stars does not significantly increase the predicted 
ionizing fluxes. The ultraviolet fluxes of OB stars can be increased by the presence of stellar winds 
(not included in the stellar models we adopted). Although winds can increase the Hell-ionizing flux by 
orders of magnitude, the corresponding increase for the H-ionizing flux in O stars has been found to be 
small [27\. The increase can be more significant for cooler B-stars |2H], but these stars do not dominate 
the ionizing photon budget in the continuous star-formation models, in which O-stars are continuously 
replenished. Furthermore, the hydrostatic stellar models typically reproduce the properties of Galactic 
HII regions |16j . so that an increase by the required factor of ~ 2 would make it more difficult to reconcile 
the models with these observations. 

Figure 2: Colour distributions of the Lya blobs. The different lines correspond to stellar population 
models with continuous star formation at an age of 10 8 yr for Z = O.2Z (dashed), Z = 0.01Z Q (dotted) 
(both computed using the models from [2^]) and Z = (solid, assuming an IMF and mass range 100 — 
500 M Q as in our analysis of LBG spectra). We have included the effect of the opacity of the intervening 
intergalactic medium. The shaded ranges indicate the range of expected colors as various properties of 
the metal-free population are varied, such as their age (between 10 5 — 10 8 yrs), the IMF (for the three 
models given in ref. |17|). and star-formation history (between a burst and continuous star formation at 
a constant rate) . All three colour histograms suggest that about 20% of the blobs need to be made of 
purely popIII stars. Note that no correction was made in the LAB colours for Lya line emission or for 
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possible continuum contribution from foreground and background objects. 
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Figure 1: 
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